Exact analytical soliton solutions in dipolar BEC 
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The bright, dark and grey solitons are well-known soliton solutions of the Gross-Pitaevskii equation for 
the attractive and repulsive BEC. We consider solitons in the dipolar BEC of the fully polarized particles, 
speaking of the dipolar BEC we mean both the magnetized BEC and the electrically polarized BEC. We 
show that these two types of the dipolar BEC reveal different behavior of the collective excitations. This 
is related to the fact that the electric and the magnetic fields satisfy to the different pairs of the Maxwell 
equation set. Thus we consider them independently. We obtain exact analytical solutions for the bright, 
dark, and grey solitons in the magnetized (electrically polarized) BEC when they propagate parallel and 
perpendicular to an external magnetic (electric) field. Comparison of spectrum of the linear collective 
excitations for the two kinds of the dipolar BEC is presented as well. 



INTRODUCTION 

Solitons as fundamental nonlinear excitations in dipolar 
BEC have been actively studied iflh- flTll . These studies are 
based on using of the generalization of the Gross-Pitaevskii 
(GP) equation for dipolar BEC H-H. Electrically po- 
larized BEC was also considered as a dielectric medium 
1 2111 . an electric field caused by polarization was explicitly 
introduced there. It has been shown that at consideration of 
the GP equation for dipolar BEC whole expression for po- 
tential energy of dipole-dipole interaction have to be used 
explicitly [22]. It allows to introduce electric field caused 
by dipoles, which satisfy to the Maxwell equations. In this 
case the GP equation for the dipolar BEC can be presented 
in an nonintegral form, and it is coupled with the pair of 
the Maxwell equations. This was developed for the electri- 
cally polarized BEC ll22ll . In this paper we will show that 
analogous representation of the GP equation for the mag- 
netized BEC can be done. We present an explicit expres- 
sion for the potential energy of magnetic dipole interaction, 
or in other words the energy of the spin- spin interaction, 
which leads to corresponding Maxwell equations. We de- 
rive our model from the first principles (the many-particle 
Schrodinger equation) using method of the quantum hy- 
drodynamics (QHD) II23I1 - 12911. which has been used for 
studying of the dipolar BEC Il30ll- ll37ll . However, dynamics 



described in review papers 111811 . [19]). 



of dipole direction had been included in previous research, 
so it was more complicated model hiding some features of 
the model, which are necessary to understand. 

In this paper we develop reduced form of the QHD equa- 
tions for fully polarized BECs. Nevertheless, even reduced 
model reveals different evolution of electrically polarized 
and magnetized BECs. Moreover, we show that nonlin- 
ear Schrodinger equations obtained in our model from the 
reduced QHD equations differ from well-known general- 
ization of the GP equation has been used for the dipolar 
BEC description (see, for instance II 1811 . IU9lP . This leads 
to different properties of collective excitations for the two 
kinds of dipolar BEC, and to difference of they properties 
from obtained using the integral GP equation, which are 



We study dispersion of collective excitations for both 
kinds of the fully polarized dipolar BECs, as for the magne- 
tized (spinor) BEC and for the electrically polarized BEC. 
We discuss differences of dispersion dependencies and rea- 
sons laying behind these differences. The reduced QHD 
leads to a nonintegral GP equations for the dipolar BECs, 
just as the reduced QHD equations can be presented in non- 
integral form instead of integral GP (see Ref. lll8ll formula 
(4.3)) and integral Euler (see Ref. [18] formula (4.9)) equa- 
tions. We show that these nonintegral equations are pow- 
erful method for finding analytical solutions for solitons in 
dipolar BECs. We consider twelve cases allowing to find 
analytical soliton solutions. We obtain solutions for the 
bright, dark, and grey solitons propagating parallel (per- 
pendicular) to external field as for the magnetized and for 
the electrically polarized BECs. 

The paper is organized as follows. In Sec. II is dedicated 
to description of our model. We briefly show the scheme 
of derivation of our model from the first principles of mi- 
croscopic theory. We describe equations of the quantum 
hydrodynamic for spinning particles (spinor BEC) and for 
particles having electric dipole moment being in the BEC 
state. We describe differences between properties of the 
two kind of BECs. We show how these equations reduce 
for the fully polarized BECs. We present and prove correct 
form of the potential energy for interaction of the electric 
dipoles, and for the spin-spin interaction, which we use for 
derivation of the QHD equations. In Sec. Ill we calculate 
and discuss dispersion of the linear collective excitations 
in the two kind of the dipolar BEC. In Sec. IV we make a 
brief review of solitons in an uniform unpolarized BEC, 
we present solutions for the bright, dark, and grey soli- 
tons, which we use in the following sections. In Sec. V we 
describe the solitons in the magnetized BEC, we consider 
propagation of the solitons parallel and perpendicular to an 
external magnetic field. In Sec. VI we study the solitons 
in the electrically polarized BEC considering propagation 
parallel and perpendicular to an external electric field. In 
Sec. VII we recapture results of our paper. 
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MODEL 

Quantum hydrodynamic equations for magnetized BEC 

The Schrodinger equation describing dynamic of N 
spinning particles being involved in the long-range spin- 
spin interaction and in the short-range interaction, and in- 
teracting with an external magnetic field is 



h lb i -°i(ext] 



ihd t 4> s (R,t) = > 




CD 



where /x is the magnetic moment of particle, = —ihdi is 
the operator of momentum, <9, is the space derivative on co- 
ordinate of zth particle, Bi( ex t) is the external magnetic field 
acting on ith particle, Uij present the short-range interac- 
tion (SRI), the Green function of the spin-spin interaction 
(SSI) has form 



(2) 



for spin matrixes s" the commutation relations are 



The first term in the right-hand side of the Schrodinger 
equation CO presents the kinetic energy of particles, the 
second term is the potential energy of magnetic moments 
in an external magnetic field B i>ext , subindex i shows that 
field Bi. ex t acts on i-th particle meaning that in the nonuni- 
form magnetic field different field acts on different parti- 
cles. The third and fourth terms describe the interparticle 
interaction, the short-range and the spin-spin interactions, 
correspondingly. 

We consider spin-1 Bose particles, thus we present the 
explicit form of the spin matrixes sf for particles with spin 
equal to 1: 






Starting from the many-particle Schrodinger equation 
(Q]l and using the QHD method we derive a set of equa- 
tions for spinning particles in the BEC state. The set of 
equations consists of three equations, they are the conti- 
nuity equation, the Euler equation giving momentum bal- 
ance, and the magnetic moment (spin) evolution equation, 
which is a generalization of the Bloch equation, the last 
one describes dynamic of single spin in an external mag- 
netic field. 



The first and basic step at getting of the QHD equa- 
tions from the many-particle wave function is the defini- 
tion of the particle concentration n(r, t) as the quantum 
mechanical average of the concentration operator in the 
coordinate presentation of the quantum mechanics h = 



Sill ^( r ~~ so we nave 



n(r,t) = / Tp*(R,t)hip{R,t)dR, 



(3) 



where dR = Yln=i * s a infinitesimally small volume 
in 3N dimensional configurational space. Next, differenti- 
ating the particles concentration (01 with respect to time we 
get the continuity equation and an explicit form of the par- 
ticle current j = nv, which is not presented here but can be 
found in papers dedicated to details of equation derivation, 
see for instance JH, Q, JH, JH, J3J. Explicit form of 
other quantum hydrodynamic al variables is also presented 
in mentioned papers. Having explicit form of the particles 
current j and magnetization M we can derive correspond- 
ing equations of evolution, which are presented below. The 
reader, who wants to know some details of the derivation 
of the QHD equation, we suggest to take following papers. 
The QHD equations for spinning particles were firstly de- 
rived in Ref. Il23ll for the charged spinning particles. Gen- 
eral technic of the derivation of the QHD equations for the 
spinless particles can be found in Ref.s ll25ll . ll2^1 . and ll35ll . 
In Ref.s 1I25I1 and ||2^1 attention was stressed on the short- 
range interaction of neutral particles in the BEC state. 



General form of hydrodynamic equations and introduction of 
magnetic field caused by magnetic moments 



The continuity equation 

d t n + d a {nv a ) 







(4) 



shows conservation of the particle number and describes 
time evolution of the particle concentration. 
The momentum balance equation 

h 2 

mn{d t + vV)/ - — d a An 
Am 



4m \ n 



-^-d^[MW a d^(— 
Amu 2 \ \ n , 



= -gnVn+M^d a (^Bi t {r,t)+ J dr'G M {r, r')AF(r\ t) 

(5) 

appears in the rather complicated form. The first group 
of terms mn(d t + vV)v = mn^ is the convective time 
derivative of the velocity field v. Next three terms are pro- 
portional to h 2 and present the quantum Bohm potential 
for spinning particles, the first two of them are the quan- 
tum Bohm potential for spinless particles corresponding 
to the GP equation for scalar macroscopic wave function. 
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The last of them arises due to the spin. If spin of all parti- 
cles are parallel to each other and to an external magnetic 
field, we can rewrite the magnetization M as M = /ml, 
where 1 is a unit vector in the direction of magnetic mo- 
ment, which is parallel to direction of the external mag- 
netic field, fi is the magnetic moment of particle, n is the 
particle concentration. In described approximation evo- 
lution of the magnetization reduces to the particle con- 
centration evolution. So, discussed term reduces to zero 



and 



li 2 d fs (n\d a d f) \) = 0, since 1 is 



a constant. The right-hand side of equation (O describes 
the force density leading to the particle evolution. The first 
term describes the short-range interaction in the first or- 
der by the interaction radius Q, S, H. The 
next group of terms presents interaction of the magnetic 
moment with the external magnetic field and interaction 
between the magnetic moments (the spin-spin interaction), 
correspondingly. 

Magnetization evolution is described by the following 
equation 



d t M a + d p \M 



2m/j \ n , 



= fie^M p (m st (v,t) + j dr'CT s (r,r')M s (r',tfj. 

(6) 

The second term d l3 (M a v 13 ) is an analog of convective 
part of the time derivative of the magnetization, which is 
also convective part of the magnetization flow. The third 
term is the quantum part of the magnetization flow, which 
is an analog of the quantum Bohm potential in the momen- 
tum balance equation (0. Let us discuss now the right- 
hand side of magnetic moment evolution equation ©. The 
terms in the right-hand side give contribution of interaction 
in time evolution of magnetization. The first term presents 
interaction of magnetic moments with the external mag- 
netic field, and the last term describes spin-spin interac- 
tion. We admit that the short-range interaction, when in- 
teraction between particles with different spin projection is 
described by the same potential, gives no contribution in 
this equation. 

Set of the QHD equations (O-© is a closed set of equa- 
tions. However these equations contain integral terms de- 
scribing the interparticle interaction. We can see, from 
equations (0I and (O, that the integral terms are written 
along with the external magnetic field, so we can expect 
that these terms describe an internal magnetic field caused 
by the magnetic moments. Thus, we can explicitly write 
the internal magnetic field 

B a (r,t) = J dr'G aP (\ r - r' \)M (r',t), (7) 

and find that this field satisfies to the Maxwell equations 
V x B(r, t) = 4vrV x M(r, t), (8) 



VB(r,i) = 0. 



(9) 



Consequently, the first term in the right-hand side of equa- 
tion ([5) can be written as M^VB , where B is the sum of 
the external and internal magnetic fields. Analogously, the 
right-hand side of the magnetic moment evolution equa- 
tion appears as (// / 'h)e a01 B 1 . As a result we have a 
set of nonintegral equation explicitly including the mag- 
netic field caused by the magnetic moments coupled with 
the equations of field ([8]) and (0. 



Reduced set of QHD equations for the fully magnetized BEC 

For the fully magnetized BEC, with no evolution of the 
magnetic moment direction, we have no change in the con- 
tinuity equation (O, the magnetic moment evolution equa- 
tion disappears in this limit, and we write down a reduced 
Euler equation 

h 2 

mn(d t + vV)/ - — d a An 
Am 



h 2 „„ / d a n ■ d^n 
4m \ n 



-gnV a n + /md a \B, (10) 



where we have used that M^VB = fml fi VB = 
/xnV(lB). 

Moreover, for the fully magnetized BEC equations of 
field represent as 



and 



VB = 0, 



V x B = Air/jVn x 1. 



(11) 



(12) 



We can also put down corresponding NLSE, which can 
be directly derived from couple of the continuity and the 
Euler equations dU and (fTOl (see Ref. II25I1 ). 



ihdMr,t) = \^-^-A+g | $(r,t) | 2 -fjW(r, t)j $(r, t), 

(13) 

where the last term —ulB describes the spin-spin inter- 
action and interaction of the magnetic moments with the 
external magnetic field. This term appears instead of the 
integral term in the generalization of the GP equation for 
dipolar BEC (see, for instance, Ref.s [18] and [19]). An 
equation analogous to (TTSl was obtained for the electrically 
polarized BEC [22]. 
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Potential energy of spin-spin interaction 

Potential energy of the spin-spin interaction was directly 
derived from the quantum electrodynamical scattering am- 
plitude of two electrons [38]. This potential energy con- 
sists of two parts. The first one is well-known appearing as 
a fraction 



U, 



dd 



d 2 - 3(dr) 2 /r 2 



(14) 



and the second term is proportional to the Dirac delta func- 
tion. However, coefficient at the delta function for the spin- 
spin interaction differs from one for the spin-spin interac- 
tion. Since, the magnetic field caused by magnetic dipole 
should satisfy to the Maxwell equations divB = and 
curlB = AncurWl, which differ from the Maxwell equa- 
tions for the electric field, the last of these equations is 
written in the quasi static approximation and in absence 
of electric currents. Finally, we write down different forms 
of the Hamiltonian of spin-spin interaction 



2)f a fP jf 2 



8vr 



S af3 5(r) 



This is explicit form of 
which we can rewrite as 



„„ - - \AnS ap 6( rij ) + V?Vf— (15) 



= - M ? M J - 5 afj A t ] — . (16) 



Explicit form of the Hamiltonian is obtained using follow- 
ing identities 

1 _ % r a rP lr 2 4tt 



+ — 6 af, 8(r), (17) 



and 



A- = -4vr5(r). 
r 



(18) 



dipoles comes as Udd = — djEji, where dj is the elec- 
tric dipole moment of the second dipole, and Ej, is the 
electric field caused by the first dipole acting on the sec- 
ond dipole. Explicit form of the dipole-dipole interaction 
energy is Udd = — d"df Vf Vf (1/r^). It can be rewrit- 
ten via notion of the Green function G al3 of dipole-dipole 
interaction U d d = —dfdfGff, where Gff = dfdf—. 

For consideration of the neutral particles having electric 
dipole moment df in the BEC state we derive a set of QHD 
equations. We starting from the many-particle Schrodinger 
equation 

ihd t ip{R,t) = Hi>(R,t), 

where R = (t±, ...,rjv) is the set of coordinate of N parti- 
cles, with the Hamiltonian 



H 



4^ I 2m, F 



d?E" 




i:, <r;,i ;(;■:; 



(19) 



where E" ext is the electric field, m., is mass of parti- 



cles, h is the Planck constant, and G\ 



d?d?l/ nj is 



the Green functions of dipole-dipole interaction discussed 
above, fTy = J7(|rj— r^- 1) is the potential of the short-range 
interaction. Let us admit that the dipole-dipole interaction 
leads to evolution of both the positions of particles and the 
directions of electric dipole moments, but in this paper we 
neglect the evolution of dipole direction. Consequently we 
will simplify the Euler equation and get a closed set of the 
QHD equations. 
Analogously for spinning particles we find the continuity 



d*n + V(nv) 



0. 



(20) 



and the Euler 



mn 



h 2 

d a An 

4m 



Quantum hydrodynamic equations for electrically polarized 

BEC 

Let us start with discussion of potential energy of elec- 
tric dipole interactions. In previous subsection we have 
met the Hamiltonian of spin-spin interaction, which dif- 
fers from usually using form of the Hamiltonian of dipole- 
dipole interaction (fT4] |. Here we are going to show that 
the potential energy of electric dipole interaction differs 
from (fT4l and (fj~5b . Potential of the electric field caused 
by an electric dipole appears as tp = — (dV)(l/r) ll39ll . 
consequently the electric field has form of E = — Vty?= 
(dV)V(l/r) or it can be rewritten in the tensor form 
E a = d /3 V /3 V Q (l/r). Energy of interaction of two 



+—d l3 (-{d a n){d?n)) = -gnd a n 
4m \n ) 

+ P l3 d a Ei t + P d a J dr'G M {r,r')P~<{r',t) (21) 

equations, where n is the particle concentration, v is the 
velocity field, P" is the polarization of the medium, or the 
density of the electric dipole moment. The last term in the 
Euler equation (|2T| | contains the integral presenting internal 
electric field caused by the polarization of medium, this 
electric field has form 

E? nt (r,t) = J dr'G^(r,r')P^r',t). 
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We can obtain that this electric field satisfies to the 
Maxwell equations 



and 



VE(r,t) = -47rVP(r,t), 



V x E(r,t) = 0. 



(22) 



(23) 



For the fully polarized BEC of particles having electric 
dipole moment the Euler equation (l2~TT l reduces to 

mn(d t +vV)v a - ^-d a An+^-d 13 ( -(d a n)(d n 
4m 4m \ n 



-gnd a n + ndd a (\E), 



(24) 



where 1 is the unit vector in direction of polarization formed 
by external field P = dnl. Along with the Euler equation 
one of the Maxwell equations reduces too, so we have 



VE(r,t) = -47rd(lV)n(r,t), 



(25) 



instead of equation (l22l . Equation d23l has no change. 

The generalization of the GP equation, for particles hav 
ing electric dipole moment, corresponding to equation 
appears as 



ihd t $(r,t) 



2m 



A+ 5 |$(r,t)| 2 -dlE(r,t)J$(r,t), 

(26) 



where the last term is proportional to — cflE, it describes 
both the interaction among dipoles and the interaction of 
dipoles with the external electric field, this term appears 
instead of the integral term (see, for instance, Ref. lfl8ll 
formula (4.3), JH formula (15), and JH formula (7)). We 
see that equation ( |26l l has the same structure as equation 
(ITS! for spinning particles, but in equation (TTSl we have the 
magnetic moment \i instead of the electric dipole moment 
d, and the magnetic field B instead of the electric field E. 
Difference in evolution of the two kinds of BEC reveals 
in the Maxwell equations for the electric E and magnetic B 
fields £□}, <H3, (l2lT >. and d25T >. Equation (|26j was obtained 
in Ref-TH- 



LINEAR EXCITATIONS 

We find that the spectrum of collective excitations in the 
fully magnetized and in the fully electrically polarized are 
different. For the electrically polarized BEC we have 



h k gn k 4irn d k cos 9 



+ 



+ 



(27) 



4m 2 m m 
where cos# = k z /k, for details of calculation see Ref. 



1 22|], and for the magnetized BEC spectrum appears as 



2 H 2 k A gn k 2 4:7rn [j, 2 k 2 sin 2 9 

to = H . (28) 

4m z m m 



Details of the calculations for the collective excitation 
spectrum of the magnetized BEC are presented in Ap- 
pendix. 

Formulas (l27l l and (l28l l are generalization of the Bogoli- 
ubov spectrum for the fully polarized dipolar BECs, for the 
BEC of particles having electric dipole moment and for the 
BEC of spinning particles, correspondingly. The first term 
in the right-hand side of equations (l27l l and (l28l l describes 
dispersion of the free quantum particle (dispersion of the de 
Broglie wave), it appears from the quantum Bohm potential 
in the Euler equation (see equation ([TUt for the magnetized 
BEC, and equation d24l) for the electrically polarized BEC). 
The second term describes the short-range interaction and 
contains the constant of short-range interaction g, which 
does not contain contribution of the dipole-dipole interac- 
tion. It depends neither on the electric dipole moment d no 
on the magnetic moment \x. The last term in formula (lZ7l 
(in formula (l28l l) appears due to interaction of the electric 
dipoles (due to the spin-spin interaction). Comparing for- 
mulas d27l) and (l28l l we see that they differ by the last term. 
Interaction of the electric dipoles leads to d 2 k 2 , where the 
external electric field is directed along the z axes, but in- 
teraction of the magnetic moments gives —^ 2 k\, where 
^1 = + ky assuming that the external magnetic field 
directed along z axes either. This difference is caused by 
the fact that the electric and the magnetic fields satisfy to 
different pairs of the Maxwell equations. 

For comparison of the last terms in the formulas (|27| 
and d28l) we show two functions £^ = x cos 2 9 and £ M = 
—x sin 2 9 on Fig.l (a) and they modules of the Fig.l (b). 

Formulas (l27T l and (l28l l differ from results have been ob- 
tained (see for instance [18] and [19]). There two results 
we have difference with. One of them obtained for the 
spinor BEC of spin-1 atoms [40]-[49]. This difference is 
related with the fact that different constants of the short- 
range interaction were used in Refs. [l40ll . [41 1 for descrip- 
tion of interaction between the species with different spin 
projection. We suppose that the short-range interaction is 
described with one interaction constant, but we include the 
long-range spin-spin interaction. Interaction of the mag- 
netic moments with an external magnetic field is also in- 
troduced in some of these papers [43] -[48]. The magnetic 
field introduces two terms, one of which is given by the 
linear Zeeman term and the quadratic Zeeman term. These 
terms give the potential energy of particles with different 
projections of the spin being in the external magnetic field. 
So, we have considered different models for the spinor 
BEC description. 

The second result was obtained on the way of general- 
ization of the scalar GP equation, including in this equa- 
tion a term describing dipole-dipole interaction I112ll - [|l4ll . 
I Hi, (Hfl. It gives results similar to (|27j and (l28l . but not 
the same. Comparison of formula (l27l l with results of other 
authors can be also found in Ref. 122]. Let us discuss 
some reasons for this difference. Firstly, the generalized 
GP equation includes interaction of the magnetic dipoles 
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FIG. 1. (Color online) The figure illustrates the dipolar part of the 
dispersion dependence for the collective excitations in the dipolar 
BECs. We picture ^ = xcos 2 6 and £jf = —a; sin 2 9. f<j depicts 
the last term in formula ((27} describing dispersion of the electri- 
cally polarized BEC. £m illustrates the last term in formula (1281 
describing dispersion in the magnetized BEC. Figure (a) shows 
£d and £fcf , Plane on figure (a) presents unpolarized BEC. Figure 
(b) presents module of ^ and £m . 



and interaction of the electric dipoles in the same way, but 
we have shown that energies of the electric dipole interac- 
tion and the spin-spin interaction have different forms. This 
reveals in difference of the Maxwell equations for the mag- 
netic field caused by the magnetic dipoles and the electric 
field caused by the electric dipoles. It also reveals in differ- 
ent properties of the collective excitation of the magnetized 
BEC and the electrically polarized BEC, as for linear ex- 
citations, their dispersion is presented by the formulas (l27l l 
and ( |28l l, and for the nonlinear soliton excitations described 
in the following sections. 



SOLITONS IN UNPOLARIZED BEC 



Following Ref.s |51], ll52ll we present analytical solu- 
tions for the bright, dark, and grey solitons in the unpolar- 



ized BEC, which will be generalized for the dipolar BECs 
in the next sections. 

The bright and dark solitons have a trivial constant veloc- 
ity field showing propagation of perturbation of the particle 
concentration. We have increasing of the particle concen- 
tration in the bright soliton existing in the attractive BEC, 
and decreasing of the particle concentration in the dark 
soliton existing in the repulsive BEC. The grey soliton ex- 
ists in the repulsive BEC and reveals soliton structure as in 
the particle concentration and in the velocity field. 

Presenting formulas appear as solutions of the following 
set of QHD equations for nonlinear plane waves propagat- 
ing in the three dimensional unpolarized BEC, so equations 
emerge in the one dimensional form 



d t n + d x (nv) = 0, 

and 

ra i a ^ h * as , h2 a ( fan) 2 
mn[Ot+vo x )v o x n-\ o x [ 

4m Am \ n 



(29) 

-gnd x n. 
(30) 



Bright soliton 

The bright soliton solution appears for the attractive 
BEC. A three dimensional attractive BEC is unstable to the 
collapse. Therefore, for consideration of the attractive BEC 
we should consider the quasi-one-dimensional BEC. How- 
ever, to get main picture of possible soliton solution, we 
can formally consider propagation of the plane soliton (one 
dimensional propagation) in the three dimensional BEC. 

We suppose that the particle concentration n(x,t) de- 
pends on the space coordinate x and time t only in the com- 
bination £ = x — V t, where V is the velocity of the bright 
soliton propagation. Assuming that perturbations caused 
by the soliton disappears at large distances from the soliton 
we come to the following boundary conditions 



lim n(£) 



0. 



(31) 



We also assume that the velocity field is an arbitrary con- 
stant v(x, t) = Vo. 

One can obtain that the particle concentration in the 
bright soliton has form of 



n(0 



2m I E 



1 



9 



cosh/ 



(32) 



where soliton amplitude n and phase a appear as 

2m I E I 



n 



9 



and 



a 



y/2 E mi 



(33) 



(34) 
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correspondingly. In the formulas (l32l-(l34t E = c — V 2 /2 
is an analog of the kinetic energy of the flow, but E has 
minus in front of the velocity square, c is a constant l5lll . 
1 12]. E < 0, g < 0. Width of the soliton D usually 
introduced as a distance between edges of the soliton at the 
half of its height. Hence, the width D equals to H ln(y2 + 
l)/{m^2\E\). 



Dark soliton 

The dark soliton exists in the unpolarized BEC with the 
repulsive short-range interaction g > 0. It reveals as an 
area of rarefication of the particle concentration. We can 
come to this solution from the set of QHD equations ( |29l 
and (]30l or corresponding one-dimensional GP equation. 
For solution we have to assume that the particle concentra- 
tion n(x, t) depends on the coordinate x and time t in the 
following combination £ = x—V Q t, hence n(x, t) = n(£), 
where V Q is a constant velocity showing velocity of the 
soliton propagation in the positive x-direction. Assuming 
that soliton gives no perturbation in the particle concentra- 
tion at the infinite distance we write the following boundary 
conditions 



lim n(£) 

x— >±oo 



n . 



(35) 



As for the bright soliton we assume that the velocity field 
is an arbitrary constant v(x, t) = V . 
The dark soliton solution appears as 



n(0 



2 / Jmgn a 
n Q tanh I £ 



(36) 



We have deep rarefication in the center of the dark soli- 
ton since tanh(0) = 0. The particle concentration in- 
creases monotonically in both directions from the center 
of the soliton. 



Grey soliton 

The grey soliton in the unpolarized BEC, as the dark soli- 
ton, exists at the repulsive short-range interaction (g > 0). 
It reveals as an area of rarefication of the particle concen- 
tration. In contrast with the dark soliton, we have a pertur- 
bation of the velocity field along with a perturbation of the 
particle concentration. 

At consideration of the grey solitons we suppose that, as 
in previous cases, the particle concentration is a function 
of £ = x — u t, so we have n(x,t) = We also 

assume that, there are perturbations of the velocity field 
v(x, t) in vicinity of a velocity Vo, and this perturbation de- 
pends on £ either. Thus, we can write velocity field v(x,t) 
as v(x, t) = Vq + Vi (£), where Vq, as in the previous sub- 
section, is an arbitrary constant current velocity associated 
with the Madelung's fluid background motion and V\ is an 
arbitrarily large current velocity perturbation. 



The particle concentration and the velocity field satisfy 
to the following boundary conditions 



and 



lim ra(£) = n , 

x— >±oo 



lim V(0=0- 



(37) 



(38) 



It was shown that velocity of the grey soliton propagation 
should satisfy to the following conditions Vq — w^ZS < 

«o<K) + y^(seeRef. jS). 

The particle concentration in the grey soliton appears as 



n(0 = n [I- A 2 



cosh 



y/ mgn A 2 



(39) 



where 



1 



m 



gn 



(uq - Vq) 2 , 



and < A 2 < 1. One can find 

(Vo - Uq)A 2 



Vi(0 



cosb/ 



2 / ^mgn A 2 



(40) 



(41) 



A 2 



for perturbation of the velocity field. 

For the future references it is useful to introduce the 
phase of the grey soliton a g as a g = ^/mgn A^/h. Area 
of the maximum rarefication of the particle concentration 
is in the center of the grey soliton n min = n(0). If A < 1 
when n(0) = n (l — A 2 ) > 0. We see that the particle 
concentration does not necessarily equals to zero, whereas 
in the dark soliton n(0) = 0. The particle concentration 
n(£) monotonically increases in both directions from the 
center of the soliton. It is interesting to consider a limit 
case when Uq = Vq, then Vi(£) = 0, A 2 = 1, and n(£) 
(]39l simplifies to the dark soliton solution (l36l . From fig- 
ures (2)-(4) we see the form of the grey soliton. Figure (3) 
shows that the velocity field perturbation appears as an area 
of increasing (decreasing) of the velocity field at Vq > u Q 
(V < n ). 



SOLITONS IN MAGNETIZED BEC 

This and next section are dedicated to solitons in the 
dipolar BEC. Magnetized BEC is described by the set of 
equations ©, (HJ, © and ([Toll. Electrically polarized BEC 
is governed by equations (l20l l. (l23l . (l24l and (|25] |. Solv- 
ing the Maxwell equations for special cases we put these 
solutions in the Euler equation. Hence, the Euler equation 
represents in a simpler form, which is an analog of the Eu- 
ler equation for the unpolarized BEC (130K Consequently, 
we can represent and use solutions, described above, for 
description of the dipolar BEC. 
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propagation of the soliton perpendicular to external 
magnetic field 

Equation ( fTTT) gives B' x = 0, where ' is the derivative 
on x, so we have B x = 0, as at infinity (in equilibrium) x 
projection of magnetic field equals to zero. Equation ([12} 
leads to two relations B' = and B' z = 4-7T jin' ', the first of 
them gives B y = 0. Derivative of B z we put in the Euler 
equation ( fTOl) . and we find 



mn(d t +vd x )v- 



■- — o x n+- — o x 
4m 4m 



(d x n) 



n 



-gMnd x n, 
(42) 



where g M = g — 47r/i 2 . This equation corresponds to the 
spectrum of linear excitations obtained above. We see that 
contribution of the magnetization reveals as effective "at- 
tractive" changing of the short-range interaction constant. 



Bright soliton in magnetized BEC 

In the unpolarized BEC the bright soliton exists at the 
attractive short-range interaction. In considering limit the 
magnetization leads to additional attracting between par- 
ticles, so magnetization of the attractive BEC does not 
break condition of the bright soliton existence. We have 
g = - | g \, consequently g M = -(| g | +47T// 2 ). We 
see that for the system under consideration the magnetiza- 
tion leads to increasing of module of g M . Only amplitude 
of the bright soliton depends on the interaction constant 
n o ~ 1/ I 5m I ■ Hence, we find decreasing of the bright 
soliton amplitude due to the magnetization of the attractive 
BEC. 



Dark soliton in magnetized BEC 

The dark soliton exists in the BEC with the repulsive 
short-range interaction (g > 0), hence the magnetization 
leads to decreasing of the effective interaction constant 
down to g rn = g — A-K^i 2 . First of all it reveals in the widen- 
ing of the dark soliton. Using of the Feshbach resonance 
allows to change the of the interaction constant g. Thus, 
at small enough positive g. Contribution of magnetization 
can prevail. That leads to disappearing of the dark soliton, 
and on stability of the repulsive BEC to the collapse. 



Grey soliton in magnetized BEC 

Grey soliton solution in the unpolarized BEC presented 
by formulas (l39T>-d4TT) . In our case we have to replace g 
with g M . It changes the phase and "amplitude" A of the 
grey soliton. To trace consequence of changing of the in- 
teraction constant we present Fig.s (2)-(4), where n(^,g) 
and Vi(£,,g) are shown. From these figures we see that 



decreasing of the interaction constant due to magnetization 
leads to decreasing of the soliton amplitudes. We have the 
decreasing of the amplitude of the particle concentration 
and the amplitude of the velocity field. Figure (4) allows 
to compare profile of the particle concentration and the ve- 
locity field. 



propagation of the soliton parallel to external magnetic field 

Let us consider a nonlinear plane perturbation propagat- 
ing parallel to external uniform magnetic field, which is 
directed along the z axes. In this case space dependence 
of hydrodynamic variables n, v and the magnetic field B 
reduces to dependence on z. Equation VB = ([Til rep- 
resents as B' z = 0, where ' means derivative on z. We 
find that magnetic moments evolution gives no influence 
on perturbation propagating parallel to the external field in 
the magnetized BEC, since B z only gives a contribution in 
the set of QHD equations © and ( fTOl) . Using equation (fT2l 
we obtain that B x and B y equal to zero. Finally, we have 
equation 

mn{o t +vo z )v—- — o z n+- — o z \ = —gno z n, 

4m 4m \ n J 

(43) 

which mathematically has no difference with equation 
d30l) . Consequently, solutions of equation d30l) presented 
by formulas (l32l l. d36l> . (l39l l and (l4~TTl are valid for our case 
and we can conclude that at solitons propagation along 
direction of the external magnetic field, the external field 
gives no influence on properties of solitons. They depend 
on the short-range interaction only. This conclusion corre- 
sponds to the properties of linear excitations in the magne- 
tized BEC d28l) . From formula (128} we see that at propaga- 
tion parallel to the magnetic field 6 = and the last term 
equals to zero. Thus, we have the unchanged Bogoliubov 
spectrum. So, nonlinear perturbation receive no influence 
of the magnetization at propagation parallel to the external 
field, as it was for the linear excitations (see section IV). 



SOLITONS IN ELECTRICALLY POLARIZED BEC 

We have described our results for the magnetized BEC. 
Let us present solitons in the fully polarized BEC of 
molecules or atoms having electric dipole moment. As we 
will see they have some different properties. 



propagation of the soliton perpendicular to external electric 

field 



Considering perturbation of hydrodynamic variables n, 
v, and E are functions of x — Vot assuming that an external 
electric field is applied along z axes. Therefore we con- 
sider propagation of the perturbation perpendicular to the 
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FIG. 2. (Color online) The figure describes the profile of the par- 
ticle concentration N = n(£) /no in the grey soliton. We consider 
the dependence of TV = n(£)/no on the reduced interaction con- 
stant G = 4n #/(mVo 2 ) and S = mVo£/(2h) at « = Vq. 




FIG. 3. (Color online) The figure describes the profile of the 
velocity field v = Vi(£)/Vb in the grey soliton on the reduced 
interaction constant G = 4n g/(mVo 2 ) and S = mVo£/(2H) at 

u = Vq. 



external field. We find that the electric field gives no contri- 
bution in the evolution of the particle concentration n and 
the velocity field v, the last term in equation (l24t equals to 
zero, since d x E z = that follows from equation (|22] |. and 
we have no difference between considered here case and 
the unpolarized BEC. Consequently, we have that soliton 
solutions (l32l l. (l36l l. ( |39l l and (l4ll propagate in the elec- 
trically polarized BEC perpendicular to the external field 
with no influence of the polarization. 



propagation of the soliton parallel to external electric field 

We now consider propagation of a nonlinear perturbation 
parallel to the external electric field, so we suppose that n, 
v, and E are functions of £ = z — v t. In this case the last 
term in equation (l24l appears as ndd z E z = —Aird 2 nd z n 
that we find using equation d22t . Thus we can combine two 
terms in the right-hand side of equation fl24) , which present 
z projection of the force density F and write it in the form 
of F z = — {g+Aird 2 )nd z n. We have that the dipole-dipole 
interaction of the fully polarized electric dipoles leads to 
change of the interaction constant. All difference, we find 
in compare with the unpolarized case, is replacement of g 
with g d = g + And 2 . Final form of the Euler equation is 



mn(d t +vd z 



h 2 o3 h 2 . 
. —d 3 z n+ —d z 
4m 4m 



(d z n) 



n 



-g d nd z n, 
(44) 



[52J we find 
and grey soli- 



Using methods described in Ref.s [51], 
changes of properties of the bright, dark, 
tons. General form of these solitons are given by formulas 
d32l) . (1361 d39l l and (gB, but instead of g we should put g d 
containing contribution of the electric dipole moment. 



Bright soliton in electrically polarized BEC 




FIG. 4. (Color online) On this figure we present the profiles of 
the particle concentration and the velocity field on one figure to 
compare they structures. They are shown as functions of G = 

An g/(rnVQ ) and H = mVb£/(2fi) at u a = Vq. 



In contrast with the magnetization, which leads to ad- 
ditional effective attraction g M = g — 4nfj> 2 at the wave 
propagation perpendicular to the external magnetic field, 
the electric polarization gives an effective repulsion g d = 
g + And 2 revealing at wave propagation parallel to the ex- 
ternal electric field. Consequently, polarization of the at- 
tractive BEC g < leads to decreasing of the module of 
the effective interaction constant g d = — \ g \ +And 2 . For 
the bright soliton it reveals in both an increasing of the soli- 
ton amplitude (l33l l and a constriction of the area of soliton 
existence. 



Dark soliton in electrically polarized BEC 

The dark soliton can also exist in the electrically polar- 
ized BEC. At theoretical description it appears as a solution 
of equation d44l) corresponding to formula (l36l . Increasing 
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of the polarization leads to increasing of the effective inter- 
action constant g d = g + 4ixd? . Consequently, the width 
of the dark soliton (g > 0) D = h/(2^mg d n ) • In ^-j- 
decreases with the increasing of polarization. The ampli- 
tude n of dark soliton is not affected by the interaction. It 
is determined by the equilibrium particle concentration. 

Grey soliton in electrically polarized BEC 

Let us describe changes of the grey soliton properties 
in the electrically polarized BEC in comparison with the 
unpolarized BEC. Fig.s (2) and (3) allow to get the conse- 
quence of polarization of the BEC on properties of the grey 
soliton. From figures (2) and (3) we see that the electric 
polarization of the BEC leads to increasing of the ampli- 
tudes of concentration and velocity field in the grey soli- 
ton. The width of the grey soliton increases as well. It hap- 
pens due to increasing of the effective interaction constant 
9d = 9 + 4vr(i 2 . 

CONCLUSIONS 



electric dipoles are presented in the paper. We obtained 
that the generalization of the Bogoliubov spectrum are 
anisotropic. Form of anisotropy is similar for both kinds 
of dipolar BEC, but they were different. They have same 
form, but direction of axes of anisotropy differs on 90 de- 
gree. Anisotropic terms caused by dipolar effects have 
differen sign. Electric dipoles lead to effective repulsion, 
since we have g + 4ird 2 cos 2 9 instead of the interaction 
constant of the short-range interaction g. Magnetic dipoles 
give effective attraction, hence we obtained g—47r/j, 2 sin 2 9 
instead of g for the unpolarized BEC. Our model allows to 
get several exact analytical soliton solutions for the dipo- 
lar BEC. We found solutions for the bright, dark, and grey 
solitons in the electrically polarized and in the magnetized 
BECs. We obtained solutions in cases when solitons prop- 
agate parallel or perpendicular to the external field. We 
discover that at propagation parallel (perpendicular) to the 
external magnetic (electric) field solitons in the magnetized 
(electrically polarized) BEC are not affected by the magne- 
tization (electrical polarization), and they have same prop- 
erties as solitons in the unpolarized BEC. In opposite limits 
we calculated contribution of the dipolar effects in the am- 
plitude and phase of solitons. 



As a summary we recapture main results of our paper. 
We studied fundamental linear and nonlinear collective ex- 
citations both in the electrically polarized BEC and in the 
magnetized BEC. For this aim we have developed the QHD 
models for the fully polarized BEC, which are simplifica- 
tion of the QHD models for the dipolar BEC including evo- 
lution of the dipole directions we had developed in earlier 
papers. These models consist of the couple of QHD equa- 
tions (the continuity and Euler equations). This couple of 
the QHD equations can be presented as a nonintegral non- 
linear Schrodinger equation, which is the generalization of 
the well-known GP equation on dipolar BECs. The QHD 
equations or corresponding nonlinear Schrodinger equa- 
tion have to be coupled with the Maxwell equations. The 
QHD equations for the electrically polarized BEC contain 
an electric field caused by electric dipoles of the medium, 
and, correspondingly, the QHD equations describing mag- 
netized BEC contain an magnetic field caused by the mag- 
netic moments. The electric field and the magnetic field 
satisfy to different pair of the Maxwell equations, conse- 
quently, we get differences in properties of the two kind 
of the dipolar BEC. These differences reveal in properties 
of the collective excitations obtained in the paper. Gener- 
alization of the Bogoliubov spectrum on the dipolar BECs 
contains a new term caused by dipolar effects. The elec- 
tric polarization gives k 2 d 2 cos 2 9, whereas the magne- 
tization leads to — k 2 fi 2 sin 2 9. These differ from well- 
known results obtained at using of the integral general- 
ization of GP equation for the dipolar BEC. These dif- 
ferences are caused by mistreating of the potential energy 
of dipole-dipole interaction. Correct forms of the poten- 
tial energy for the spin- spin interaction and interaction of 



APPENDIX 

We describe details of the calculation of the spectrum of 
collective excitations for the fully magnetized BEC. From 
the continuity equation dU in the linear approximation, as- 
suming that n = n + 5n, v = + <5v, and B = B + SB, 
where n Q and B are equilibrium values of the particle con- 
centration and the magnetic field, equilibrium value of the 
velocity field equals to zero, 5n, 6v, and 5B are small per- 
turbations of the equilibrium state, for the Fourier ampli- 
tudes of the hydrodynamic variables, we find 

5n = n -^—. (45) 

The Euler equation ( fTQb in the linear approximation repre- 
sents as 

h 2 

— iumn 5\ + zk/c 2 - — 5n = —gn i\i5n + unoi)a6B z , 

4m 

(46) 

multiplying this equation on k and using equation (l45l l 
in the first term we get relation between the perturbation 
of particle concentration and z-projection of the magnetic 
field, which is 

h 2 

— uj 2 m5n + k 4 - — 5n = —gn k 2 5n + nn k 2 5B z , (47) 

4m 

Equation dTzb gives three scalar equations 

k y SB z — k z SB y = 4Trk y fj,5n, (48) 

k x 5B z — k z 5B x = 47rk x fi5n, (49) 
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and 



h y 5B x 



0. 



(50) 



From the equation (l50l we can get 5B X via 5B y , using it 
and algebraic form of equation ([Til , which is 



k x 5B x + kySBy + k z 5B z = 0, 
we obtain 5B y via 5-B z in the following form 



5B V 



k k 

^5B Z . 

k 2 ± 



(51) 



(52) 



Next, we put (1521 in equation (1481 1 and find relation for 
8B Z , which is a part of the Euler equation (l46l . via the 
perturbation of the particle concentration as follows 



5B Z 



k 2 



Sn, 



(53) 



where k\ 



kl + ky, and k 2 



kl + k 2 y + k 2 z . There- 



fore, we can write k\ = k 2 sin 2 9, where is the angle 
between direction of wave propagation given by k and di- 
rection of the external magnetic field B . Putting ( 1531 ) in 
the linearized Euler equation ( |47l we get formula i 



t 
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